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Motivated by observation of very high tunnel magnetoresistance (TMR) in Fe-MgO-Fe magnetic 
tunnel junction devices, we propose a theoretical model for these devices based on a single-band tight- 
binding approximation. An effort is made to capture the band dispersions over the two dimensional 
transverse Brillouin zone. In the transport direction, spin dependent Hamiltonian is prescribed for 
Ai and A5 bands. Non-equilibrium Green's function formalism is then used to calculate transport. 
Features like voltage dependence of TMR are captured quantitatively within this simple model and 
the trends match well with the ones predicted by ab-initio methods and experiments. 

PACS numbers: 72.25.-b, 85.75.-d, 75.47.-m, 75.47.Jn, 85.35.-p 



Magnetic tunnel junction (MTJ) devices have emerged 
as one of the candidates for random access memory appli- 
cations. The prediction of high tunnel magnetoresistance 
(TMR) for crystalline MgO barrier P, [3| was followed by 
observations of upto 200% TMR ratios at room temper- 
ature in Fe-MgO-Fc and CoFc-MgO-CoFe MTJ devices 
[1, Since then, there has been an increased effort to 
integrate them into practical devices. Although, ab-initio 
P, M, H, 01 and empirical tight binding studies have 
been reported, their computational complexity limit their 
use for rapid device prototyping. 

In Rcf. 0], we report parameters for bcc Fc(lOO), 
which reflect the band dispersions for majority and mi- 
nority spin bands over the two dimensional transverse 
Brillouin zone within a simple single-band tight-binding 
(SBTB) method Q. The transmission calculated using 
this method matches well with the one calculated us- 
ing EHT-NEGF method (extended Hiickel theory, non- 
equilibrium Green's function). In this paper, we propose 
parameters for MgO as shown in Table I and complement 
them with the Fe parameters in Ref. ^ to capture the 
band structure effects in Fe-MgO-Fe MTJ devices. This 
approach provides a simple and computationally inex- 
pensive platform to explore MTJ devices. Furthermore, 
it gives an inherently simple and intuitive understanding 
for the underlying device physics. 

The schematic of a MTJ device and the band dia- 
gram is shown in Fig. 1(a) and Fig. 1(b) for parallel 
(P) and anti-parallel (AP) configurations respectively. In 
P configuration, the magnetizations of the two contacts 
are in the same directions and in AP configuration, the 
magnetizations are opposite. Therefore, the P config- 
uration has higher current densities than AP configura- 
tion. This change in current densities defines the TMR as 
(Jp — Jap)/ Jap- Higher TMR signifies a higher signal- 
to-noise ratio and hence is desirable. In this paper, we 
normalize TMR as (Jp — Jap)/{Jp + Jap) to compare 
our results with Refs. [1, [1] . 



For quantum transport, we use NEGF formalism [T^j- 
For each band, we start with a one dimensional single- 
band tight-binding Hamiltonian 0] for the device shown 
in Fig. 2: 



H: 
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where Ei,o is the band offset and Ul is the Laplace po- 
tential linearly dropped across the insulator region. For 
MgO, the hopping parameter to is given in Table I and 
to for Fe bands are given in Ref. Q . Each device lattice 
point in Fig. 2(c) corresponds to a unit cell in Fig. 2(b) 
[shown by the dotted box] . The resulting dispersion is of 
the form e{k) = E^o + 2to[l — cos(fca/)], where a; = 4.2A 
is the cubic lattice constant for MgO. At Fe-MgO inter- 
face, on-site elements are taken as average of to for Fe 
and MgO. The off-diagonal elements arc taken such that 
the resulting Hamiltonian is Hcrmitian as needed to en- 
sure that the energy eigenvalues are real and current is 
conserved. The Green's function is then calculated as: 



G^[{Ei+l{)+)I -HsBTB-^c\ 



where Ec = —toe^'^"'""'^ and c=[l,2] for the left and 
right contact respectively. The transmission is given 
as T{Ei) = tr{riGT2G^), which is fc|| independent and 
transmission per unit area is calculated as: 
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where a = 2.86A is the Fe cubic lattice constant. Current 
density in P and AP configuration for each spin orienta- 
tion is then computed as: 



J = j dEi TsBTB [fi — 72] 



(4) 
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(a) Parallel Configuration (P) 




(b) Antiparallel Configuration (AP) 
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FIG. 1: A magnetic tunnel junction (MTJ) device. An insulating tunneling barrier is sandwiched between two metallic 
ferromagnetic layers, (a) Parallel (P) configuration along with its energy band diagram showing a cosine dispersion obtained 
using single-band tight-binding method. The magnetization of the two contacts are in the same direction, (b) Anti parallel 
(AP) configuration is shown, where the magnetizations of the two contacts are in opposite directions. 




FIG. 2: The device structure for transport calculations, (a) 
A schematic of the device used, (b) The device structure 
used in ab-initio calculations, reproduced from Ref. [1]. The 
unit cell is shown by dotted line which corresponds to one 
lattice point in SBTB-NEGF method, (c) Real space lattice 
for SBTB-NEGF calculations. First lattice points in Fe are 
treated as interface for which onsite Hamiltonian element is 
sum of hopping parameters of Fe and MgO and off diagonal 
elements are those of Fe or MgO depending on the neighboring 
lattice point. 



TABLE I: Parameters for MgO. Ut is the MgO barrier height, 
to is the hopping parameter. The parameters for Fe are in Ref. 
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FIG. 3: Current density for Ai and A5 bands in Fe-MgO-Fe 
MTJ device. Current densities for parallel and anti-parallel 
configurations using SBTB-NEGF method. The parameters 
for MgO are shown in Table I and for Fe are in Ref. [^. 



For Ai band, Ut is taken to be half of the band gap in 
local density approximation (LDA) of the density func- 
tional theory [ll|. Although LDA underestimates band 
gap, we use this value to compare our results with Refs. 
[a, |6|. Ub for A5 band is used as a fitting parameter, to 
for both the bands are estimated by fitting current levels 
to Ref. @. The calculated current densities for Ai and 
A5 bands using these parameters are shown in Fig. 3. 

In Fig. 4, the total P (Jp-totai) and AP (JAP-totai) 
current density is shown which match well with ab-initio 
results 0. Jp-totai is dominated by Ai band current 
density. This is due to the lower potential barrier seen 
by the Ai band. On the other hand, the total AP cur- 
rent density is dominated by A5 band till about IV. After 
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FIG. 4: Total current densities for P and AP configurations. 
The total current densities using SBTB-NEGF method for P 
and AP configurations are shown by solid and dotted lines 
respectively. The crosses are data from Ref. . The low bias 
transport in P configuration is dominated by Ai band and in 
AP configuration is dominated by As band. The high bias 
transport in both configurations is dominated by Ai band. 



this voltage, there is a sharp increase in AP current of 
Ai band and this band current density starts to domi- 
nate the total current density. This is due to half-metallic 
nature of Ai band. To elaborate more, any current in a 
particular configuration is made up of two components - 
the majority spin current J-\ and the minority spin cur- 
rent J I . In anti-parallel configuration, both these current 
channels are made up of a minority band in one contact 
and a majority band in the other as shown in the band 
diagram in Fig. 1(b). At voltage below about IV, cur- 
rent is small because the band edge of one of the contacts 
is still not in the bias window. Once the applied bias is 
high enough, the band is pulled within the bias window 
and the current starts to increase. In Fig. 5, TMR calcu- 
lated using SBTB-NEGF method is shown. TMR values 
match well with those obtained using the ab-initio model 
in Ref. [1]. The rapid decrease in TMR after about IV 
is also consistent with other ab-initio studies for a three 
layer device 0. The TMR calculated using SBTB-NEGF 
method also follows a similar bias dependent trend as re- 



ported in experiments [l^]. To the best of our knowledge, 
such a bias dependence and quantitative agreement has 
not yet been captured within a simple single-band tight- 
binding model or an effective mass based model. 

We have presented a single-band tight-binding model 
for bee Fe-MgO-Fe magnetic tunnel junctions in [100] 
direction. We have tried to capture the band structure 
effects by using band parameters for different symmetry 
bands in the contacts and the barrier region. Features in 
TMR which are manifestation of the electronic structure 
of material were captured quantitatively within a simple 
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FIG. 5: Tunnel magnetoresistance (TMR). TMR calculated 
using SBTB-NEGF method is shown, x's show TMR calcu- 
lated using ab-imtio method reproduced from Ref. The 
bias dependence of TMR is captured well within this simple 
single-band tight-binding model. 



model. 
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